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1.0  INTRODUCTION 

The  contribution  of  skin  friction  drag  to  total  drag  is  frequently  required  in 
transonic  aerodynamic  testing.  This  occurs,  for  example,  when  total  drag  measurements 
are  made  in  tunnel  tests  where  flight  Reynolds  numbers  cannot  be  duplicated.  In  such 
cases,  the  friction  drag  corresponding  to  test  Reynolds  number  must  be  extracted  from 
the  tunnel  data  and  replaced  by  the  friction  drag  corresponding  to  flight  Reynolds 
number.  Such  data  manipulations  must  be  accomplished  by  the  application  of  correlation 
or  prediction  techniques.  Therefore,  the  objective  of  the  work  described  herein  was  to 
determine  a suitable  theoretical  procedure  for  the  calculation  of  skin  friction  drag  in 
transonic  turbulent  flows  at  arbitrary  Reynolds  numbers.  Of  primary  interest  is  the 
accuracy  of  tire  calculation  method;  however,  it  would  be  advantageous  if  the 
computation  technique  were  also  fast  and  required  a minimal  amount  of  computational 
resources. 

Boundary  layers  over  two-dimensional  planar  and  axisymmetric  bodies  in  transonic, 
turbulent  flow  were  computed  using  five  theoretical  calculation  techniques  and  the  results 
are  compared  with  experimental  data.  In  addition  to  skin  friction,  other  boundary-layer 
parameters  are  compared  to  provide  a more  complete  test  of  the  computational 
techniques. 

The  calculation  methods  chosen  for  evaluation  are  presented  in  Table  l in  order  of 
increasing  complexity  with  respect  to  basic  formulation  and  turbulence  modeling.  Also 
shown  in  Table  1 are  the  machine  requirements  (in  terms  of  the  AEDC  IBM  370/165) 
and  extended  capabilities  (second-order  boundary-layer  effects)  pertaining  to  each 
calculation  procedure.  The  capabilities  of  each  method  are  described  in  more  detail  in 
Section  2.0. 

Experimental  results  were  chosen  from  two  investigations:  a planar,  supercritical 
airfoil  experimental  study  (Ref.  1),  and  an  investigation  involving  an  axisymmetric, 
waisted  body  of  revolution  (Ref.  2).  Both  experimental  studies  were  concerned  with  the 
measurement  of  mean  How  turbulent  boundary-layer  properties  at  transonic  speeds  and 
adiabatic  wall  conditions. 


2.0  CALCULATION  METHODS 

The  calculation  techniques  evaluated  can  be  classed  in  three  categories  which  are  the 
integral  boundary-layer  equations,  finite  difference  boundary-layer  equations,  and  the 
ensemble  averaged,  compressible,  Navicr-Stokes  equations.  The  computational  techniques 
used  are  documented  elsewhere  and  only  brief  descriptions  of  the  equations  solved  and 
the  numerics  and  inputs  required  arc  given. 
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Table  1. 


Requirements  end  Extended  Capabilities  of  Computational  Methods 


c 

Computer  Requirements 

Extended  Capabilities 

X 

Memory 
Required , 

Typical 

CPU  Tima, 

Relative 

Preparation 

Time 

r 

Method 

^ame 

Method 
Type  ^ 

K-bytas , 
<IBMJ70/165 

Min, 

IBM  370/165 

Prep.  Time 
for  Method  1 

Turbulence 

Model 

Transverse 

Curvature 

Longitudinal 

Curvature 

Kree-Stremn 

Turbulence 

Vail 

Roughness 

Wall 

Transpiration 

Trans i tion 

Wliir  field 

Integral 

Boundary 

Layer 

40 

1 

Algebraic 

Herring- 

ttellor 

Finite 

Difference 

Boundary 

Layer 

130 

0,5 

1.5 

Algebraic 

X 

X 

X 

X1 

AALK  M 

Finite 

Difference 

Boundary 

Layer 

165 

6.0 

1.5 

V -Equation 
of  Turbulence 

X 

X 

X2 

tfi  lenx 

Finite 

Difference 

Boundary 

Layer 

265 

7.0 

3 

2-Equatlnn 
ot  Turbulenra 

X 

X 

X 

X 

X3 

Del  vert 

Finite 

Difference 

Navier- 

Stokea 

460 

100 

20 

Algebraic 

X 

X 

X 

X1 

1 Artificial— user  can  specify  the  relative  proportions  of  laminar  and  turbulent  viscosities  with  axial  variation 
*"  Location  of  tramfitlcnt  can  be  controlled  by  specifying  amount  of  TKE  in  the  layer 
^ Same  as  2 with  addition  of  controlling  transit inn  hy  also  specifying  turbulence  length  scale 


AEDC-TR-79-12 


A6DC-TR-79-12 


2.1  INTEGRAL  BOUNDARY-LAYER  METHOD,  WHITFIELD  PROGRAM 


The  integral  boundary-layer  technique  used  is  described  in  Ref.  3.  Not  unlike  other 
integral  methods,  this  technique  involves  transforming  the  partial  differential  equations  of 
the  boundary  layer  into  ordinary  differential  equations  in  terms  of  integral  parameters, 
i.e.,  the  displacement  thickness,  5*.  momentum  thickness,  d,  etc.  The  approach  is  to 
simultaneously  solve  the  momentum  integral  equation 


d . _ 5*  dUe 

- *•'’«>  * r — • 
dx 


Ue  dx 


(i) 


and  the  mean-flow  kinetic  energy  integral  equation 

1 a , |3 


2rk  o U3  dx 

W r c C 


d i.  S JU 

— -r 

U dx 


where 


D 


(2) 


(3) 


(4) 


oc 


Equations  (1)  and  (2)  are  valid  for  two-dimensional  planar  (k  = 0)flow,  and  axisymmetric 
(k  = 1 ) flow  with  the  terms  involving  the  transverse  curvature  effects  neglected.  The  main 
advantage  of  this  method  over  otheT  integral  techniques  is  attributed  to  two  recent 
developments:  (1)  the  improved  anlaytical  description  of  the  compressible,  turbulent 
boundary-layer  velocity  profile  presented  in  Ref.  4,  where  it  is  shown  that  the  entire 
turbulent  boundary  layer  can  be  reasonably  well  described  by  one  analytical  expression 
involving  the  shape  factor,  H,  skin  friction  coefficient,  Cf  e,  and  Reynolds  number  based 
on  boundary-layer  momentum  thickness,  Refl,  and  (2)  the  improved  velocity-temperature 
expression  for  turbulent  boundary  layers  with  nonunity  Prandtl  numbers,  developed  in 
Ref.  5. 


To  describe  the  effect  of  turbulence  within  the  boundary  layer  (i.e.,  to  model  the 
Reynolds  stress  term),  a three-layer,  algebraic,  eddy  viscosity  model  is  used.  In  the  region 
nearest  the  wail  (y+  < 100),  the  total  shear  stress  (molecular  plus  turbulent)  is  assumed 
constant  and  equal  to  the  wall  value.  The  eddy  viscosity  in  the  middle  region  (y+  > 100, 
Cm  < e o ) is  given  by  (Ref.  6) 
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€ 


m 


(0-4-ly) 2 


(?u 

dy 


In  the  outer  region  (defined  for  y greater  than  the  minimum  y where  e0 
constant  eddy  viscosity  as  suggested  by  Clauser  (Ref.  7)  is  used: 


(6) 

em),  a 


£„  = 0.0168  Ue5* 


(7) 


where 


The  integral  equations  of  the  boundary  layeT  are  then  cast  into  finite  difference 
form  and  solved  by  the  predictor-corrector  method  used  by  Nash  (Ref.  8).  Inputs 
required  for  the  program  include  initial  values  of  shape  factor,  momentum  thickness, 
free-stream  Mach  and  Reynolds  numbers,  and  the  axial  variation  of  edge  Mach  number 
and  body  surface  radius  (for  axisymmetric  flow).  For  the  comparisons  presented  herein, 
the  measured  values  of  shape  factor  and  momentum  thickness  at  the  first  measurement 
station  were  used  as  initial  conditions  from  which  the  initial  value  of  Cf>e  is  determined 
from  a correlation  of  CfiE  = Cf,e(Refl,H)  (Ref.  3). 

2.2  FINITE  DIFFERENCE  BOUNDARY-LAYER  METHODS 


2.2.1  Harring-Mellor  Program 


The  method  of  Herring  and  Mellor  (Refs.  9 and  10)  encompasses  the  solution  of  the 
compressible,  two-dimensional  (planar  or  axisymmetric),  laminar  and/or  turbulent 
boundary-layer  equations.  The  equations  solved  are 


<9  d 

-(rpu)  + ^-{rpV)  = 0 


(9) 


du  d u d 1 1 e 1 d 

57 + pv  w = U«-T,“  * 7 57  (rro)  ( 1 0) 


„ l?ll„  , <?h„ 

Pu  ~ + Pv  — 

px  ay 


r r)  y 


(ID 


which  are  valid  for  either  laminar  or  turbulent  flow  if  r0  and  q0  are  defined  as  the  total 
shear  stress  and  total  heat  flux,  respectively, 


ro/P  * 


<11  ' V ' > 


(12) 
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tJli 


qo/P  = k^--<v'h'> 


(13) 


The  gross  effects  of  turbulence  are  simulated  by  using  the  "effective  viscosity  hypothesis" 
of  Mellor  (Ref.  11).  The  effective  viscosity  for  the  total  boundary  layer  is  given  by 


c>(.X  R) 

".i-M 


and  the  effective  conductivity  by 


where 


a...  = U.  S 


V V 

1 / vef  \ 

R<TP  ,ipf,\u/  H«gv 

IT 

o.4iv  p 


sv 


SU. 


R = 


(14) 


(15) 


(16) 


(17) 


</>  = 1 + -, r 

X3  + C6.9)3 


x = RX 


(18) 


® = X (X  <0.016),  = 0.016  (X> 0.016) 


(19) 


Considerable  flexibility  is  possible  in  terms  of  the  specification  of  boundary 
conditions.  Provisions  have  been  made  to  include  the  effects  of  wall  transpiration  and 
roughness  as  well  as  the  capacity  of  treating  an  adiabatic  wall  and/or  variable  wall 
temperature  case.  An  artificial  method  of  simulating  transition  can  be  effected  by 
specifying  the  amount  and  location  of  turbulent  viscosity  to  be  linearly  combined  with 
the  molecular  viscosity  (Ref.  10). 

To  effect  a solution,  the  boundary-layer  equations  (Eqs.  9 through  11)  are  first 
transformed  into  an  intermediate  form  by  using  a modified  version  of  the,  Probstein-Elliot 
transformation  (Ref.  12).  Then,  by  utilizing  axial  direction  finite  differences,  the 
equations  are  transformed  again  into  ordinary  differential  equations.  The  equations  are 
then  written  for  each  grid  point  across  the  boundary  layer  and  the  resulting  characteristic 
matrix  is  solved  by  Gaussian  elimination.  Initial  velocity  and  temperature  profiles  (either 
laminar  or  turbulent),  boundary-layer-edge  velocity  distribution,  and  any  wall  conditions 
pertaining  to  the  specification  of  wall  temperature  distribution  (or  an  adiabatic  wall),  wall 
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roughness,  or  transpiration  are  required  for  initialization  of  the  program.  (It  is  worthy  to 
note  that  the  incompressible  version  of  the  program  received  a "good"  rating  from  the 
1968  "Boundary  Layer  Olympics11  held  at  Stanford  University  (Ref.  13.)).  In  this  case, 
the  initial  turbulent  velocity  profile  was  generated  by  using  the  results  of  Ref.  4 and  the 
first  measured  values  of  Cf>e,  H,  and  d. 

2.2.2  ABLE-II  Program* 


Similar  to  the  equations  solved  by  the  Herring-Mellor  method,  the  ABLE-II  program 
("ABLE"  referring  to  "Algorithm  for  Boundary  Layer  Type  Equations”  of  Patankar  and 
Spalding  (Ref.  14))  obtains  a finite-difference  solution  to  the  compressible, 
two-dimensional  (planar  or  axisymmetric),  laminar  and/or  turbulent  boundary-layer 
equations,  taking  into  account  any  effects  of  wall  transpiration  and  transverse  curvature 
(longitudinal  curvature  and  wall  roughness  effects  are  neglected).  Calculations  can  also  be 
made  with  an  adiabatic  wall  and/or  variable  wall  temperature.  In  terms  of  nondimensional 
variables,  the  equations  solved  by  the  ABLE-II  method  are 

d 3 

— (p  r 50  — tp  rv)  = 0 (20) 

ax  ay 


r = r w + y cos  0 


(24) 


‘This  program  was  developed  by  Dr.  i.  A.  Bcnek,  ARO,  Inc.,  and  is  noi  reported  elsewhere. 
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The  plus  sign  in  Eq.  (24)  is  for  external  flow,  the  minus  sign  for  internal  flow. 


Computation  of  the  Reynolds  stress  goes  beyond  the  simple  algebraic  eddy  viscosity 
concept  in  that  the  following  particular  partial  differential  equation  (Ref.  15)  describing 
the  turbulent  kinetic  energy  (TKE)  across  the  boundary  layer 


where 


(25) 


(26) 


is  solved  in  conjunction  with  the  boundary-layer  equations,  Eqs.  (20  through  22). 
Laminar  flow  can  be  represented  by  setting  the  TKE  across  the  boundary  layer  equal  to 
zero.  Transition  to  turbulent  flow  can  be  simulated  by  initially  setting  the  TKE  to  some 
specified  low  value  and  allowing  the  turbulence  to  grow.  It  should  be  pointed  out  that 
although  reasonably  good  mean-flow  quantities  can  be  computed  through  transition,  the 
location  of  transition  is  not  known  a priori.  The  transition  location  can  be  adjusted  by 
changing  the  input  value  of  the  maximum  TKE  in  the  layer. 

The  numerics  of  the  program  involve  solving  the  boundary-layer  equations  in 
normalized  von  Mises  coordinates  by  means  of  the  Patankar-Spalding  finite  difference 
procedure  (Ref.  14).  Inputs  required  to  obtain  a solution  are  free-stream  Mach  and 
Reynolds  numbers,  axial  pressure  and  wall  temperature  distributions,  wall  transpiration 
rates,  and  the  initial  starting  velocity  and  temperature  profiles,  the  latter  being  either 
laminar  or  turbulent.  For  the  comparisons  presented  herein,  the  ABLE-1I  solutions  were 
begun  from  near  the  stagnation  point  and  the  transition  location  was  varied  such  that  the 
calculated  and  measured  values  of  momentum  thickness  at  the  first  measuring  station 
were  approximately  equal. 

2.2.3  Wilcox  Program 

The  Wilcox  program  (Ref.  16)  solves  the  compressible,  two-dimensional  (planar  or 
axisymmetric)  boundary-layer  equations  for  laminar  and/or  turbulent  flow.  The  equations 
can  be  written  as 

d , , id 
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where  q is  given  either  by 


q = 


f 

pT 


or  by  simultaneously  solving  Eqs.  (27  through  29)  with  the  equation 


dq  4 (9  32  1 d 

<" 5* pv T * ~ 5 T ~ eT ’ T T 


(28) 


(29) 


(30) 


The  eddy  diffusivity,  e,  in  Eqs.  (29  through  30),  is  assumed  to  be  related  to  the  TKE  by 


1 (31) 

<0 


where  cj  is  the  turbulent  dissipation  rate. 

The  turbulence  model  utilized  is  more  sophisticated  than  those  described  previously 
in  that  two  equations  are  solved  for  the  turbulent  viscosity,  the  TKE  equation  plus  an 
equation  describing  the  turbulent  dissipation  rate  (TDR).  The  equations  used  by  Ref.  16 
for  describing  the  TKE  and  TDR  are 
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where  0,  0*,  a,  and  cr*  are  constants,  $ = e1/2  co,  and  a and  a*  are  functions  of 
molecular  viscosity,  TK'E,  and  TDR. 

Flows  with  boundary  conditions  such  as  wall  roughness,  transpiration,  varying  wall 
temperature  and/or  an  adiabatic  wall  can  be  computed  with  relative  ease. 

In  order  to  solve  Eqs.  (27  through  33),  the  physical  boundary-layer  coordinates  are 
transformed  to  the  computational  plane  by  the  Probstein-Elliot  (Ref.  12)  and  Levy-Lees 
(Ref.  17)  transformations.  The  resulting  block  tridiagonal  system  of  algebraic  equations 
is  then  solved  by  the  method  of  Blottner,  et  al.  (Ref.  18). 

Inputs  required  for  the  program  are  the  usual  external  conditions  such  as  free-stream 
Reynolds  number,  variations  of  the  boundary-layer-edge  Mach  number,  axial  wall 
temperature  distribution  (or  an  adiabatic  wall),  etc.,  plus  a constant  value  of  the 
free-stream  turbulent  kinetic  energy  and/or  turbulence  length  scale.  From  the  latter 
stipulation,  it  is  necessary  to  begin  all  calculations  from  an  initial  laminar  starting  profile. 
Similar  to  the  ABLE-1I  program,  transition  to  turbulent  flow  can  be  controlled  by 
varying,  the  input  values  of  free-stream  TKE  and/or  the  free-stream  turbulence  length 
scale.  For  the  comparisons  presented  herein,  the  Wilcox  solutions  were  begun  from  near 
the  stagnation  point  and  the  transition  location  was  varied  such  that  the  calculated  and 
measured  values  of  momentum  thickness  at  the  first  measuring  station  were 
approximately  equal. 

2.3  FINITE  DIFFERENCE  NAVIER-STOKES-DEIWERT  PROGRAM 

The  Deiwert  program  (Refs.  19  through  21)  solves  the  time-dependent, 
ensemble-averaged  (Refs.  22  and  23)  compressible,  two-dimensional,  Navier-Stokes 
equations  for  arbitrary,  two-dimensional  planar  geometries.  It  should  be  noted  that  since 
the  equations  have  been  ensemble-averaged,  time-dependence  corresponds  here  to  a period 
much  longer  than  that  associated  with  turbulent  fluctuations,  e.g.,  see  Ref.  22.  The 
program  is  basically  an  extended  and  improved  version  of  the  MacCormack  method  (Ref. 
24)  originally  introduced  in  1969  (Ref.  25).  The  present  program,  however,  casts  the 
Navier-Stokes  equations  into  integral  form,  whereas  MacCormack 's  formulation  is  in 
differential  form  (Ref.  26).  Written  in  vector-integral  form,  Deiwert  formulates  the 
time-dependent,  compressible,  mean-flow,  Navier-Stokes  equations  as 


^ /ud(vol),/  H 


ndg 


0 


(34) 
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where 


V = 


uex  + vev 


H = 
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lu  dv\ 
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br  + <9x  / 

(40) 

and  e£,  e£  are  unit  vectors  in  the  x and  y directions,  respectively,  and  if  is  a unit  vector 
normal  to  a grid  cell  face. 


A simple,  two-layer,  eddy  viscosity  turbulence  model  (Ref.  22)  is  employed  within 
the  boundary-layer  region  to  simulate  the  effects  of  turbulent  fluctuations  on  the  mean 
flow.  For  boundary-layer  flow,  the  turbulent  viscosity  is  assumed  to  be 


€ 


av 

t?x 


(41) 


where 


r 


l = O.iy 


— exp 


for  the  inner  region  and 


26 


Jv- 


f = o.oys 


(42) 


(43) 


for  the  outer  region. 
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In  principle,  any  physically  plausible  boundary  conditions  such  as  adiabatic  or 
variable  temperature  walls  and  wall  transpiration  can  be  implemented  (except  for  subscale 
wall  conditions  such  as  roughness).  To  simulate  transition  from  laminar  to  turbulent  flow, 
the  total  viscosity  (moleeufar  plus  turbulent)  is  set  equal  to  the  molecular  viscosity 
upstream  of  a specified  axial  location  on  the  body.  Downstream  of  this  location,  the 
eddy  viscosity  (computed  from  the  algebraic  model)  is  added  to  the  molecular  viscosity 
in  the  normal  fashion. 

The  numerical  method  casts  the  equations  into  an  integral,  primitive  variable, 
conservative  form  (e.g.,  see  Ref.  27)  and  utilizes  the  predictor-corrector,  finite  difference 
method  of  MacCormack  (Refs.  20  and  24)  for  solution.  The  method  is  impulsively 
started  from  uniform  free-stream  conditions  imposed  throughout  the  computational 
domain.  The  solution  then  evolves  with  time  until  a steady  state  is  reached. 

Because  u rather  extensive  nonorthogonal  grid  (Ref.  20)  must  be  generated  before 
the  program  can  be  implemented,  a considerable  amount  of  effort  is  required  to  ensure 
that  important  features  of  the  flow  will  be  properly  resolved.  Thus,  in  general,  the  inputs 
to  the  program  are  significantly  more  cumbersome  than  those  needed  for  previously 
described  boundary-layer  methods.  For  example,  in  order  to  generate  a mesh  around  an 
airfoil,  one  must  specify  where  and  how  the  grid  cells  arc  distributed  around  the  leading 
edge  and  along  the  chord.  In  addition,  to  gain  computational  efficiency  in  the  normal 
direction  to  the  surface  (Ref.  28)  the  mesh  is  split  into  two  or  more  regions.  In  regions 
where  large  flow  gradients  are  anticipated,  the  grids  describing  these  regions  are 
exponentially  stretched  in  order  that  acceptable  numerical  resolution  can  be  achieved.  In 
any  case,  creation  of  the  proper  grid  is  not  a trivial  matter.  Therefore,  care  must  be 
exercised  in  allocating  the  number  and  location  of  mesh  points,  the  amount  of  stretching, 
etc.,  such  that  all  pertinent  features  of  the  flow  are  reasonably  well  resolved  numerically. 

3.0  EXPERIMENTAL  DATA 

A literature  search  revealed  that  very  little  reliable,  experimental,  transonic, 
two-dimensional,  turbulent  boundary-layer  data  have  been  published.  Only  two 
appropriate  experimental  investigations  were  found.  Both  experimental  studies  were 
conducted  in  the  Royal  Aircraft  Establishment's  (RAE)  transonic  facilities. 

3.1  RAE  AIRFOIL  (PLANAR) 

The  experimental  study  of  Cook,  McDonald,  and  Firmin  (Ref.  1)  was  conducted  in 
the  RAE  8-ft  by  6-ft  (2.44  m by  1.83  m)  slotted  wall,  transonic,  continuous  flow  wind 
tunnel  and  involved  the  probing  of  the  turbulent  boundary  layer  over  a planar, 
supercritical  airfoil  (chord  = 0.61  m.  maximum  thickness  = 73.76  mm),  shown  in  Fig.  1, 


15 


AEDC-TR-7B-12 


at  high  subsonic  free-stream  Mach  numbers  and  at  various  angles  of  attack  and  Reynolds 
numbers.  The  flow  conditions  ranged  from  subcritical  (Miocai  < 1)  at  negative  angle  of 
attack  to  supercritical  (Mioca|  > 1)  at  positive  angle  of  attack.  Data  acquired  were 
surface  pressure  distributions  and  mean-flow  boundary-layer  profiles  deduced  from 
traverses  of  pitot  and  static  pressure  probes.  Oil  flow  was  also  utilized  to  detect  the 
location  of  boundary-layer  separation.  In  the  cases  to  be  used  herein,,  boundary-layer 
transition  was  artificially  induced  on  the  upper  surface  at  three  and  eleven  percent  chord 
by  surface  roughness.  It  was  determined  that  at  some  high  Mach  and  Reynolds  numbers 
the  roughness  height  (762  (i  at  the  high  Mach  number  cases,  127  ju  for  the  lower  Mach 
numbers)  was  larger  than  would  normally  be  used.  The  tests  were  conducted  such  that 
temperature  equilibrium  was  achieved  between  the  model  and  flow  (adiabatic  wall).  For 
the  cases  to  be  presented,  the  tunnel  pressure  fluctuations  were  reported  to  be  less  than 
1.2  percent  of  the  free-stream  dynamic  pressure.  Skin  friction  was  inferred  from  the 
velocity  profile  measurements  by  utilizing  the  equivalent,  incompressible  logarithmic  form 
of  the  velocity  profile  given  by  Winter  and  Gaudet  (Ref.  29).  Boundary-layer  integral 
properties  were  defined  such  that  any  effects  of  pressure  gradient  normal  to  the  surface 
could  be  accounted  for  (e.g.,  see  Ref.  30). 


Surface  pressures  were  the  only  measured  quantities  for  which  uncertainties  were 
determined.  At  the  low  Reynolds  number  (2.7  x 106),  the  uncertainty  in  Cp  was 
reported  to  be  less  than  ±0.0064;  at  the  high  Reynolds  number  (6.5  x 106),  less  than 
±0.0026.  The  airfoil  angle  of  attack  and  tunnel  Mach  number  settings  were  within  ±0.01 
deg  and  0.001,  respectively.  Wall  interference  corrections  were  not  included  in  the 
presentation  of  the  data. 

3.2  WAISTED  BODY  OF  REVOLUTION 

The  wind  tunnel  investigation  conducted  by  Winter,  Rotta,  and  Smith  (Ref.  2)  using 
an  axisymmetric,  waisted  body  of  revolution  involved  the  detailed  probing  of  the 
turbulent  boundary  layer  as  it  developed  over  the  body  at  high  subsonic  to  low 
supersonic  free-stream  Mach  numbers.  The  main  objective  of  their  study  was  to  produce 
an  axisymmetric.  converging  flow  with  an  adverse  pressure  gradient  in  order  to  investigate 
the  effects  of  Mach  number,  pressure  gradient,  and  streamline  convergence  and  divergence 
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on  the  development  of  the  turbulent  boundary  layer.  The  model,  shown  in  Fig.  2,  was 
1.524  m in  length  and  had  a maximum  diameter  of  0.254  m.  The  boundary  layer  was 
tripped  near  the  nose  of  the  body  at  approximately  x/L  = 0.025.  The  tests  were 
conducted  such  that  adiabatic  wall  conditions  were  achieved.  Surface  pressure 
measurements,  mean-flow  boundary-layer  velocity  profiles  deduced  from  pitot  pressure 
traverses,  and  surface  skin  friction  determined  by  the  razor  blade  technique  (Ref.  31) 
were  obtained.  Measurements  of  tunnel  free-stream  turbulence  were  not  reported. 
Boundary-layer  integral  properties,  defined  by 


=/K)0-£> 


(44) 


(45) 


were  evaluated  by  assuming  the  static  pressure  was  constant  through  the  boundary  layer. 
The  body  radius  term  (1  + y/rw)  was  retained  in  the  definitions  so  that  the  momentum 
integral  equation  could  be  used  as  a check  on  the  validity  of  the  experimentally 
determined  boundary-layer  integral  parameters. 


It  was  noted  in  Ref.  2 that  the  validity  of  surface  skin  frictions  obtained  by  the 
razor  blade  technique  is  questionable  since  it  has  not  been  fully  validated  in  transonic 
flows  with  pressure  gradients  (also,  see  the  comments  of  Rubesin,  et  al.,  in  the 
introduction  of  Ref.  32). 

Surface  pressure  coefficients  were  reported  to  be  within  ±0.005  for  a Reynolds 
number  of  107  while  skin  friction  coefficients  were  estimated  to  be  accurate  to  within 
±0.0002,  assuming  that  the  calibration  procedure  (Ref.  31)  is  valid.  The  model  shape  was 
correct  to  within  0.01  in.  (0.254  mm)  of  the  design  shape  while  the  slope  was 
constructed  to  within  0.001  in./in.  (0.0254  mm/ mm)  of  design. 
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4.0  RESULTS  AND  DISCUSSION 


4.1  RAE  AIRFOIL 

The  experiments  of  Ref.  1 with  the  RAE  2822  supercritical  airfoil  (Fig.  1)  represent 
commonly  encountered  two-dimensional,  transonic,  airfoil-type  flow  Fields  (both 
subcritical  and  supercritical).  Boundary-layer  measurements  were  only  made  on  the  upper 
surface.  Three  sets  of  measurements  were  considered,  denoted  as  Cases  2,  9,  and  12.  Case 
2 (M^  = 0.676,  Re^,.  = 5.7  x 1 06 , a = -2.18  deg)  represents  a fully  subcritical  flow  with 
generally  mild  favorable  and  adverse  pressure  gradients  on  the  upper  surface.  Case  9 (M„ 
= 0.73.  Re.^  = 6.5  x 106-.  a = 3.19  deg)  and  Case  12  (M^  = 0.73.  Re^,.  = 2.7  x 106  , a 
= 3.19  deg)  are  representative  supercritical  airfoil  flows  with  severe  adverse  pressure 
gradients  caused  by  shocks  on  the  airfoil's  upper  surface. 

4.1.1  Pressure  Distributions 

The  measured  airfoil  surface  pressure  distributions  are  shown  in  Figs.  3a,  b,  and  c for 
Cases  2,  9,  and  1 2.  respectively.  Also  shown  are  the  surface  pressures  computed  by  the 
Deiwert  program  (recall  that  the  external  flow  field  is  calculated  by  the  Nuvier-Stokes 
method,  whereas  the  boundary-layer  methods  require  information  about  the  external  flow 
field  as  input). 

Figure  3a  illustrates  that  agreement  between  computed  and  measured  surface 
pressures  for  Case  2 is  not  good,  particularly  on  the  upper  surface  where  the  computed 
results  tend  to  exhibit  oscillations  upstream  of  about  50-percent  chord  possibly  caused  by 
relatively  large  axial  grid  spacing.  Lower  surface-  calculations  are  somewhat  smoother  but 
fall  short  of  accurately  computing  the  extent  of  expansion  on  the  initial  portion  of  the 
airfoil.  In  contrast,  Fig.  3b  shows  good  agreement  between  computed  and  measured 
airfoil  surface  pressure  distributions.  The  calculated  shock  location  and  strength  agree 
remarkably  well  with  the  measured  data.  Lower  surface  calculations  again  fall  below  the 
quantitative  expansion  of  the  measured  data.  The  measured  and  computed  airfoil  surface 
pressures  at  lower  Reynolds  number,  shown  in  Fig.  3c.  illustrate  fair  agreement  between 
experiment  and  calculation.  On  the  upper  surface,  the  shock  location  was  computed  to 
be  downstream  of  the  experimental  location.  Also,  the  expansion  peak  indicated  by  the 
experimental  data  at  approximately  seven-percent  chord  is  underpredicted  by  the 
calculations  as  are  the  lower  surface  computations  over  the  entire  chord  length.  Good 
agreement  is  achieved,  however,  between  the  calculations  and  experiments  on  the  upper 
surface  downstream  of  65-perccnt'  chord. 
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^e„,c  = 2.7  x 10®,  a = 3.19  (case  12) 
Figure  3.  Concluded. 
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4.1.2  Boundary-Layer  Characteristics,  Case  2 

Results  from  the  five  calculation  methods  are  compared  with  experimental  skin 
friction  data  for  Case  2 in  Fig.  4a.  The  four  boundary-layer  calculation  methods  generally 
compare  well. with  the  experimental  data.  Surface  skin  frictions  as  computed  by  Deiwert's 
method  are  considerably  higher  than  those  measured,  particularly  on  the  downstream 
portion  of  the  airfoil  surface  where  the  pressure  gradient  becomes  adverse  (Fig.  3a). 
However,  it  should  be  pointed  out  that  the  number  of  grid  points  within  the  boundary 
layer  was  at  least  one  order  of  magnitude  lower  than  those  used  by  the  boundary-layer 
methods  which  could  account  for  the  observed  discrepancies  in  skin  friction  as  calculated 
by  the  Deiwert  program. 


o Experimental  Data 
Herring-Mellor 


x/c 

a.  Skin  friction 

Figure  4.  Upper  surface  boundary-layer  properties  of  the  RAE  2822  airfoil 
for  = 0.676,  ReM(C  = 5.7  x 106,  a = -2.18  deg  (case  2). 

Figure  4b  shows  the  measured  and  computed  comparisons  of  upper  surface 
boundary-layer  displacement  thickness  (5*).  Excellent  agreement  between  experiment  and 
all  calculation  methods  exist  in  the  upstream  region  where  the  pressure  gradient  is 
favorable.  In  the  adverse  pressure  gradient  region  the  displacement  thickness  increases 
rapdily  and  slightly  more  accurate  values  of  5*  are  calculated  by  the  Herring-Mellor 
method  than  any  of  the  other  boundary-layer  methods.  At  about  70-percent  chord, 
results  as  given  by  the  Wilcox  method  start  to  deviate  from  the  other  computations  and 
continue  to  calculate  values  of  5*  below  the  data  for  the  remainder  of  the  airfoil  length. 
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As  was  noted  with  the  airfoil  surface  pressure  distribution  calculated  by  the  Deiwert 
program  (Fig.  3a),  the  displacement  thickness  computed  by  this  method  exhibits 
oscillations  downstream  of  50-percent  chord.  It  is  fortuitous  that  the  experimental  data 
and  the  Deiwert  computed  value  of  displacement  thickness  agree  at  95-percent  chord. 


o Experimental  Data 


x/c 

b.  Displacement  thickness 
Figure  4.  Continued. 

Figure  4c  illustrates  better  overall  agreement  between  calculated  and  measured 
momentum  thicknesses  than  was  obtained  for  displacement  thicknesses  in  Fig.  4b.  Similar 
to  the  displacement  thickness,  the  momentum  thickness  is  seen  to  slowly  increase  in  the 
favorable  protion  of  the  surface  pressure  gradient  and  then  increases  rapidly  over  the 
remaining  50  percent  of  the  airfoil  surface  in  the  adverse  pressure  gradient  region.  The 
Herring-Mellor  and  ABLE-II  methods  yield  good  agreement  with  the  data  along  the  entire 
chord  length,  whereas  results  as  computed  by  Whitfield  and  Wilcox  methods  are  slightly 
under  the  data  aft  of  approximately  70-percent  chord.  Again,  results  calculated  by  the 
Deiwert  program  contain  oscillations  in  the  solution  and  a significant  increase  in  0 occurs 
past  90-percent  chord. 

Shown  in  Fig.  4d  is  the  airfoil  upper  surface  variation  of  the  measured  and 
computed  shape  factor.  The  data  indicate  almost  constant  values  of  this  parameter  over 
the  first  80  percent  of  the  airfoil  chord  as  do  the  calculations  by  the  four  boundary-layer 
methods.  Quantitatively,  however,  the  Whitfield  method  shows  the  overall  best  agreement 
with  the  experimental  data  which  may  be  a consequence  of  using  the  value  of  H at  the 
first  measuring  station  as  input.  As  should  be  expected  from  Figs.  4b  and  c,  the  shape 
factor  distribution  as  computed  by  the  Deiwert  method  shows  oscillations  not  found  in 
the  data  or  the  other  calculated  results. 
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c.  Momentum  thickness 
Figure  4.  Continued. 
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d.  Shape  factor 
Figure  4.  Concluded. 
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4.1.3  Boundary- Layer  Characteristics,  Case  9 

The  upper  surface  streamwise  variation  of  measured  and  computed  skin  friction  is 
shown  in  Fig.  5a.  Upstream  of  the  50-percent  chord  position  the  skin  friction  decreases 
slightly  with  increasing  x/c.  Because  of  the  presence  of  a shock  located  at  about 
55-percent  chord  and  the  accompanying  severe  compression,  the . experimentally 
determined  skin  friction  drops  significantly.  All  the  calculation  methods  exhibit  good 
qualitative  trends  with  respect  to  the  data.  Immediately  downstream  of  the  shock  at 
approximately  60-percent  chord,  all  calculation  methods  show  a decrease  and  then  an 
increase  in  skin  friction  with  the  ABLE-II,  Wilcox,  and  Deiwert  methods  showing  the 
most  severe  variations.  There  is  a slight  trailing-edge  separated  flow  region  computed  by 
the  Deiwert  technique  which  is  not  evident  from  the  experimental  data.  Also,  the 
ABLE-II  method  predicted  separation  at  about  82-percent  chord  which  was  the  only 
boundary-layer  method  to  indicate  such  a result. 


a.  Skin  friction 

Figure  5.  Upper  surface  boundary- layer  properties  of  the  RAE  2822  airfoil 
for  M„  = 0.73,  Re„,iC  = 6.5  x 106,  a = 3.19  deg  (case  9). 
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The  experimental  and  calculated  variation  of  displacement  thickness  on  the  airfoil 
upper  surface  at  the  Case  9 flow  conditions  are  shown  in  Fig.  5b.  Upstream  of  the  shock 
location  (x/c  = 0.55),  the  boundary-layer  methods  do  a good  job  of  calculating  the  values 
of  displacement  thickness,  blit  the  Deiwcrt  method  gives  results  that  are  too  high.  On  the 
downstream  portion  of  the  surface,  the  presence  of  the  shock  and  the  associated  severe 
adverse  pressure  gradient  cause  a thickening  and  retardation  of  the  boundary  layer,  thus 
causing  the  abrupt  increase  in  the  displacement  thickness.  Agreement  between  the 
experiment  and  theory  is  not  good  at  x/c  > 0.6. 


b.  Displacement  thickness 
Figure  5.  Continued. 

Figure  5c  presents  the  airfoil  upper  surface  streamwise  variation  of  momentum 
thickness.  Upstream  of  the  shock,  the  momentum  thickness  exhibits  a gradual  increase. 
At  and  downstream  of  the  shock,  a sharp  increase  in  0 occurs:  Results  from  the 
boundary-layer  methods  are  in  excellent  agreement  with  the  data  upstream  of  the  shock 
location,  whereas  the  momentum  thicknesses  as  computed  by  the  Deiwert  method  in  this 
region  arc  too  high  compared  with  the  data.  Downstream  of  the  shock,  the  abrupt 
increase  in  experimentally  determined  0's  is  computed  to  be  in  dose  agreement  to  the 
d3ta  by  the  ABLE-II  program;  however,  this  method  predicted  separation  at  x/c  = 0.8 
which  was  not  experimentally  measured  for  this  case.  The  other  three  boundary-layer 
methods  yield  results  which  lie  significantly  below  the  experimental  data  downstream  of 
the  shock,  whereas  the  Deiwert  method  calculates  momentum  thicknesses  which  arc  too 
high  but  with  the  same  rate  of  increase  as  the  experimental  data. 
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c.  Momentum  thickness 
Figure  5.  Continued. 


Shown  in  Fig.  5d  are  the  airfoil  upper  surface  shape  factor  distributions.  An  abrupt 
increase  in  shape  factor  occurs  in  the  shock  region  followed  by  a decrease  in  shape 
factor.  This  trend  is  qualitatively  computed  by  all  the  calculation  methods  except  the 
ABLE-1I  method.  The  Whitfield  and  Wilcox  methods  are  in  reasonably  good  agreement 
with  one  another  along  the  airfoil  and  are  in  good  agreement  with  experimental  data 
upstream  of  the  shock.  The  agreement  between  calculated  and  measured  shape  factors 
downstream  of  the  shock  is  not  good. 

o Experimental  Dala 

Herrinq-Mellor 

Wilcox 


Whitlield 

ABLE  II 


d.  Shape  factor 
Figure  5.  Concluded. 
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4.1.4  Boundary- Layer  Characteristics,  Case  12 

For  unknown  reasons,  the  Wilcox  method  could  not  be  made  to  run  Case  1 2.  Case 
12  is  similar  to  Case  9,  except  the  pressure  gradient  is  less  severe  and  the  Reynolds 
number  is  lower. 

Figure  6a  presents  the  measured  and  calculated  skin  friction  on  the  RAE  airfoil 
upper  surface  for  Case  12.  Unlike  the  higher  Reynolds  number  case  (Case  9),  a more 
gradual  decrease  in  surface  skin  friction  was  measured  in  the  shock  region,  0.4  < x/c  < 
0.7.  As  with  the  previous  two  flows,  no  separation  was  indicated  experimentally. 
Upstream  of  x/c  = 0.45,  the  Whitfield  method  is  in  good  agreement  with  the  measured 
skin  friction  distribution,  falls  below  that  measured  in  the  region  0.5  < x/c  < 0.7,  and 
recovers  to  fair  agreement  over  the  downstream  30  percent  of  the  airfoil  surface.  The 
Herring-Mellor  method  calculates  values  of  skin  friction  higher  than  the  measured  data 
upstream  of  the  shock,  but  is  in  good  agreement  with  measured  data  in  and  downstream 
of  the  shock  region.  The  ABLE-I1  method  gives  good  qualitative  agreement  with 
measured  data,  but  computes  skin  frictions  which  are  higher  than  the  measured  data 
downstream  of  the  shock.  It  is  interesting  to  note  that  the  ABLE-II  method  predicted 
separation  at  the  higher  Reynolds  number  Case  9 but  not  at  the  lower  Reynolds  number 
associated  with  Case  12  indicating  a possible  anomaly  concerning  this  method,  or  the  use 
of  it.  Values  of  skin  friction  calculated  by  the  Deiwert  program  are  too  high  upstream  of 
the  shock,  again  possibly  caused  by  an  inadequate  number  of  grid  points  within  the 
boundary-layer  region.  Also,  because  the  shock  location  was  predicted  downstream  of 
that  indicated  by  the  data  (see  Fig.  3c),  the  calculated  drop  in  skin  friction  associated 
with  the  severe  adverse  pressure  gradient  is  also  moved  downstream  relative  to  the 
boundary-layer  calculations.  It  is  also  noted  from  Fig.  6a  that  a small  separated  region 
was  computed  by  the  Deiwert  method  immediately  behind  the  shock  and  also  very  near 
the  airfoil  trailing  edge. 

Figure  6b  presents  the  streamwise  variation  of  displacement  thickness  over  the  upper 
airfoil  surface.  Similar  to  the  high  Reynolds  number  case  (Case  9),  a slow  increase  in 
displacement  thickness  occurs  upstream  of  the  shock.  Because  of  the  strong  adverse 
pressure  gradient  caused  by  the  shock,  a substantial  increase  in  8*  occurs  over  the  aft 
portion  of  the  airfoil.  The  boundary-layer  methods  calculate  5*  reasonably  well  except 
over  the  aft  portion  of  the  airfoil,  x/c  > 0.7.  The  Deiwert  program  does  not  compute  the 
distribution  of  <5*  well  for  this  case.  However,  it  should  be  emphasized  that  the  Deiwert 
program  is  concerned  with  computing  the  entire  flow  field,  whereas  the  other  methods 
treat  only  the  boundary-layer  region.  The  Deiwert  program  would  require  prohibitive 
amounts  of  machine  storage  and  computation  times  in  order  to  have  the  same  numerical 
resolution  (particularly  within  the  boundary  layer)  as  the  boundary-layer  methods. 
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x/c 

a.  Skin  friction 

Figure  6.  Upper  surface  boundary-layer  properties  of  the  RAE  2822  airfoil 
for  = 0.73,  Re^c  = 2.7  x 10®,  a = 3.19  deg  {case  12). 

Illustrated  in  Fig.  6c  are  the  measured  and  computed  values  of  momentum  thickness 
on  the  airfoil  upper  surface  at  the  Case  12  flow  conditions.  All  of  the  boundary-layer 
calculation  methods  give  results  which  agree  reasonably  well  with  measured  data.  Si  nilar 
to  the  results  for  Cases  2 and  9,  these  results  indicate  that  the  boundary-layer  methods 
do  a better  job  of  predicting  momentum  thickness  than  they  do  displacement  thickness. 

Figure  6d  shows  the  experimental  and  computed  airfoil  upper  surface  shape  factor 
distributions  for  the  Case  12  flow  conditions.  Unfortunately,  there  are  insufficient 
experimental  data  to  resolve  the  region  where  the  predictions  indicate  a peak  in  H.  The 
Whitfield  method  gives  good  agreement  upstream  of  the  shock  primarily  because  H is  an 
initial  condition.  Downstream  of  the  shock  the  ABLE-II  method  gives  the  best  qualitative 
agreement  with  measured  data,  but  the  quantitative  agreement  is  not  particularly  good. 
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b.  Displacement  thickness 
Figure  6.  Continued. 


x/c 


c.  Momentum  thickness 
Figure  6.  Continued. 
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o Experimental  Data 


d.  Shape  factor 
Figure  6.  Concluded. 

Shown  in  Table  2 are  comparisons  of  calculated  results  on  the  RAE  airfoil  in  the 
form  of  qualitative  judgements  as  to  wliich  calculation  method  gave  the  most  accurate 
predictions  for  each  boundary-layer  property  examined.  The  "X's"  represent  a choice  of 
the  most  accurate  method  while  the  "O's"  signify  the  method  is  slightly  less  accurate 
than  the  most  accurate  method. 

Table  2 illustrates  that  the  Whitfield  method  gave  the  most  accurate  skin  frictions 
for  Cases  9 and  12  (supercritical).  The  Herring-Mellor  program  predicted  the  most 
accurate  skin  frictions  for  Case  2 (subcritical).  Overall,  equal  accuracy  was  achieved 
between  the  Whitfield  and  Herring-Mellor  methods.  However,  the  fact  that  both  the 
Whitfield  and  Herring-Mellor  solutions  were  started  from  initial  conditions  at  the  first 
measuring  station,  and  the  ABLE-1I  and  Wilcox  solutions  were  started  upstream  of  this 
point  with  initial  conditions  such  that  the  momentum  thickness  of  the  first  measuring 
station  was  approximately  matched,  could  influence  how  well  these  methods  agreed  with 
experimental  data  downstream  of  this  point.  Unfortunately,  solutions  could  not  be 
obtained  from  all  methods  which  satisfied  the  same  set  of  initial  conditions,  which  of 
course,  would  be  the  basis  for  proper  comparison. 
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Table  2.  Relative  Accuracy  of  the  Computational  Methods  for  the 
RAE  Airfoil  Comparisons 
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Legend:  X Judged  to  be  most  accurate  method 

0 Judged  to  be  slightly  less  accurate  then  above 

- Judged  to  be  significantly  lesa  accurate  than  either  of  the  above 
* Wilcox  method  could  not  be  made  to  compute  this  flow 

4.2  WAISTED  BODY  OF  REVOLUTION 

The  waisted  body  of  revolution  model  used  for  the  experiments  reported  in  Ref.  2 
was  designed  such  that  an  axisymmetric.  converging  flow  with  an  adverse  pressure 
gradient  could  be  produced  and  the  effect  of  Mach  number,  pressure  gradient,  and 
streamline  convergence  and  divergence  on  the  turbulent  boundary  layer  could  be 
investigated.  Two  sets  of  experimental  data  taken  on  the  waisted  body  of  revolution  were 
selected  for  comparisons  with  the  boundary-layer  solutions:  a subsonic  case  with  M,,  = 
0.597,  Re„,,L  = 9.98  x 106,  and  a supersonic  case  with  = 1.398.  Re„>L  = 10.08  x 
106.  The  Deiwert  method  used  could  compute  only  planar  flows. 

Figure  7 shows  the  surface  pressure  distributions  at  the  two  selected  flow 
conditions.  The  distribution  of  pressures  upstream  of  the  body  waist  (x/L  = 0.7)  are 
similar  for  the  two  Mach  numbers  considered  in  that  strong  favorable  and  adverse 
pressure  gradients  are  produced  in  both  cases.  However,  downstream  of  the  waist  the 
subsonic  flow  tends  to  rc-expand,  whereas  the  supersonic  case  initially  begins  to  expand 
and  then  recompresses  over  the  Anal  protion  of  the  body. 
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Sym  ^od  ^eoo,  L x ^ 
o 0.597  9.98 

• 1.398  10.08 


Figure  7.  Pressure  distributions  on  the  waisted  body  of  revolution. 

4.2.1  Subsonic 

Illustrated  in  Fig.  8a  are  the  measured  and  computed  body  surface  skin  friction 
distributions  for  the  case  in  Fig.  7 denoted  by  the  open  symbols.  In  the  region  0.3  < x/L 
^ 0.7,  the  adverse  pressure  gradient  coupled  with  the  flow  convergence  associated  with 
the  thinning  of  the  body  produces  a generally  decreasing  CfiC  with  the  minimum  value 
occurring  near  the  body  waist.  The  favorable  pressure  gradient  aft  of  x/L  = 0.7  promotes 
an  increase  in  the  measured  skin  friction.  All  four  calculation  methods  give  results  that 
compare  favorably  with  the  experimental  data  in  terms  of  qualitative  trends,  although  all 
methods  quantitatively  compute  values  of  skin  friction  higher  than  the  data.  The 
Whitfield  method  yields  the  best  agreement  with  the  experimental  data.  The  values  of 
skin  friction  computed  by  the  Wilcox  program  are  high  over  most  of  the  body  but  drop 
to  near  that  measured  in  the  waist  region.  In  the  region  0.4  < x/L  < 0.55,  the  values  of 
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skin  friction  as  calculated  by  the  Hcrring-Mellor  program  are  increasing,  whereas  the  other 
calculations  (and  also  the  experimental  data)  show  the  opposite  trend.  It  is  interesting  to 
note  that  although  the  Whitfield  method  yields  the  best  agreement  with  these 
experimental  data,  the  transverse  curvature  terms  in  the  equations  are  neglected  in  this 
analysis.  However,  the  maximum  value  of  8*/rv.  is  about  0.16  and  hence  the  effect  of 
transverse  curvature  should  not  be  too  significant  for  these  flow  conditions.  The 
sensitivity  of  the  Whitfield  method  to  initial  conditions  is  illustrated  in  Fig.  21  of  Ref.  3 
for  these  experimental  data.  Four  sets  of  initial  conditions  were  used  in  Ref.  3.  All  four 
solutions  rapidly  converge,  indicating  little  sensitivity  to  initial  conditions  for  these  flow 
conditions. 

o Experimental  Data 


a.  Skin  friction 

Figure  8.  Boundary-layer  properties  on  the  waisted  body  of  revolution 
for  = 0.597,  Re^,L  = 9.98  x 106. 

Figure  8b  presents  the  distribution  of  displacement  thickness  variations  at  the 
subsonic  flow  conditions.  In  the  region  0.4  < x/L  < 0.7,  the  flow  convergence  and 
adverse  pressure  gradient  promote  a rapid  thickening  of  the  boundary  layer  and  an 
accompanying  increase  in  the  displacement  thickness.  Over  the  last  30  percent  of  the 
body,  the  flow  divergence  along  with  the  favorable  pressure  gradient  produce  a thinning 
of  the  boundary  layer  along  with  a decrease  in  the  displacement  thickness.  Generally, 
good  agreement  between  measured  and  calculated  values  of  5*  is  obtained  over  the  entire 
body  with  the  Whitfield  method  possibly  yielding  the  best  results.  The  maximum  value  of 
8*  located  at  the  waist  is  calculated  well  by  the  Herring-Mellor  method  which  computes 
values  of  displacement  thickness  higher  than  the  data  elsewhere. 
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o Experimental  Data 

—  Herring-Meilor 

—  Wilcox 

Whitfield 


b.  Displacement  thickness 
Figure  8.  Continued. 

The  measured  and  computed  axial  variations  of  momentum  thickness  are  shown  in 
Fig.  8c.  The  experimental  data  exhibit  the  same  increasing-decreasing  trend  as  the 
displacement  thickness  distribution  (Fig.  8b)  which  corresponds  to  the  successive  flow 
convergence-divergence  and  the  adverse-favorable  pressure  gradient  effects.  The 
Herring-Meilor  program  calculates  values  of  momentum  thickness  significantly  higher  than 
the  measured1  values  of  momentum  thickness  over  the  entire  body  length.  For  this 
particular  case,  results  computed  by  the  ABLE-1I  and  Whitfield  methods  show  good 
agreement  with  the  experimental  data.  Values  of  6 as  computed  by  the  Wilcox  program 
are  generally  lower  than  the  experimental  data. 

Figure  8d  gives  the  calculated  and  experimentally  determined  streamwise  variation  of 
shape  factor  at  the  subsonic  flow  conditions.  Values  of  boundary-layer  shape  factor 
calculated  by  the  Whitfield  program  show  fairly  good  agreement  with  the  data,  whereas 
the  remaining  three  methods  yield  a wide  disparity  of  shape  factor  distributions  as 
compared  to  the  experimental  data  (recall  that  H and  0 are  inputs  to  the  Whitfield 
program). 
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o Experimental  Data 
Herrlng-Mellor 


c.  Momentum  thickness 
Figure  8.  Continued. 


d.  Shape  .factor 
Figure  8.  Concluded. 
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4.2.2  Supersonic 


The  streamwise  variation  of  surface  skin  friction  for  the  supersonic  flow  conditions 
as  measured  and  calculated  is  shown  in  Fig.  9a.  Similar  to  the  subsonic  case,  flow 
convergence  along  with  the  adverse  pressure  gradient  in  the  region  0.4  < x/L  < 0.7 
produce  a gradual  decrease  in  measured  values  of  surface  skin  friction;  the  minimum 
value  again  occurs  near  the  body  waist.  Downstream  of  the  waist,  the  flow  diverges  as  the 
body  diameter  increases  and  there  is  an  increase  in  the  measured  values  of  surface  skin 
friction.  Upstream  of  the  body  waist  (x/L  = 0.7),  the  computations  are  generally  too  high 
compared  to  that  which  was  measured.  However,  values  of  Cf>e  computed. by  the 
Whitfield  method  are  in  better  agreement  with  the  data  than  are  the  other  three 
calculation  techniques.  Downstream  of  the  body  waist,  the  Whitfield  method  computed 
values  of  Cfie  which  eventually  fall  below  those  measured,  whereas  the  other  three 
methods  all  yield  skin  frictions  which  arc  higher  than  those  experimentally  determined 
except  near  the  end  of  the  body.  Over  the  last  ten  percent  of  the  body,  however,  the 
Herring-Mellor,  ABLF-fL  and  Wilcox  programs  all  calculate  skin  frictions  close  to  the 
experimental  data,  whereas  those  computed  by  the  Whitfield  method  are  low.  However,  a 
compatibility  check  performed  by  Winter,  Rotta,  and  Smith  (Ref.  2),  which  consisted  of 
using  measured  values  of  Cfe  and  H in  the  momentum  equation,  led  them  to  the 
conclusion  that  the  skin  friction  was  too  large  for  x/L  > 0.7.  They  attributed  this  to  the 
technique  used  to  measure  skin  friction  (Ref.  31). 
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a.  Skin  friction 

Figure  9.  Boundary-layer  properties  on  the  waisted  body  of  revolution 
for  = 1 .398,  Re„,L  = 10.08  x 10*. 
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Figure  9b  shows  the  measured  and  calculated  displacement  thickness  distribution 
over  the  body  at  the  supersonic  flow  conditions.  The  thickening  of  the  boundary  layer 
enhanced  by  flow  convergence  and  an  adverse  pressure  gradient  upstream  of  the  body 
waist  reflect  a corresponding  increase  in  the  measured  values  of  5*,  the  maximum  of 
which  occurs  again  near  the  waist  of  the  body.  Downstream  of  this  location,  divergence 
of  the  flow  thins  the  boundary  layer.  General  trends  are  predicted  well  by  all  four 
methods,  but  upstream  of  the  body  waist,  values  of  6*  calculated  by  Whitfield  and 
Herring-Mellor  methods  are  in  better  agreement  with  the  data  than  either  of  the  other 
methods.  Downstream  of  the  body  waist,  all  four  calculation  methods  compute  values  of 
5*  fygher  than  the  data. 


b.  Displacement  thickness 
Figure  9.  Continued. 

Shown  in  Fig.  9c  are  the  values  of  momentum  thickness  at  the  supersonic  flow 
conditions.  Upstream  of  x/L  = 0.7  (body  waist),  the  Herring-Mellor  method  tends  to 
compute  values  of  0 which  are  too  high,  whereas  those  calculated  by  the  Wilcox  program 
are  too  Low.  In  this  same  region,  the  Whitfield  and  ABLE-II  methods  calculate  values  of 
momentum  thickness  which  are  in  better  agreement  with  the  experimental  data. 
Downstream  of  the  body  waist,  the  Herring-Mellor,  ABLE-II,  and  Whitfield  methods  give 
values  of  0 higher  than  those  actually  measured.  In  this  region,  the  Wilcox  program  more 
accurately  computes  the  correct  quantitative  values  of  the  momentum  thickness. 

The  axial  distribution  of  boundary-layer  shape  factor  as  measured  and  computed  at 
the  supersonic  flow  conditions  are  shown  in  Fig.  9d.  In  this  case,  the  Whitfield  and 
Wilcox  methods  give  the  best  descriptions  of  the  shape  factor  distribution  upstream  of 
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the  body  waist  (x/L  = 0.7).  Aft  of  x/L  = 0.7,  shape  factors  computed  by  the  Whitfield 
and  Wilcox  methods  are  too  high,  whereas  those  calculated  by  the  Herring-Mellor  and 
Able-II  methods  are  in  better  agreement  with  the  experimental  data. 


c.  Momentum  thickness 
Figure  9.  Continued. 


d.  Shape  factor 
Figure  9.  Concluded. 
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' Table  3 presents  qualitative  results  as  to  which  technique  was  judged  to  yield  the 
most  accurate  boundary-layer  calculations  using  the  experimental  data  as  the  standard  of 
accuracy.  In  both  the  subsonic  and  supersonic  cases,  the  Whitfield  method  gave  the  most 
accurate  computations  of  skin  friction.  For  these  cases,  the  Whitfield  program  also  more 
accurately  calculated  almost  all  of  the  boundary-layer  integral  properties.  However,  the 
same  caution  offered  in  Section  4.1.4  concerning  initial  conditions  applies  here. 
Therefore,  the  conclusions  are  applicable  only  to  the  particular  experimental  data  used 
for  comparison,  and  the  manner  in  which  initial  conditions  were  implemented  for  each 
calculation  procedure. 

Table  3.  Relative  Accuracy  of  the  Computational  Methods  for  the 
Waisted  Body  of  Revolution  Comparisons 
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Legend:  X Judged  to  he  the  most  accurate  method 

0 Judged  to  be  slightly  less  accurate  than  above 

- Judged  to  be  significantly  less  accurate  than 
either  of  the  above 

5.0  CONCLUSIONS 

Numerical  calculations  have  been  made  to  determine  the  accuracy  of  five  computer 
programs  to  compute  the  skin  friction  on  a two-dimensional  wing  and  a body  of 
revolution  in  transonic  flow.  It  was  shown  that  of  the  five  computational  methods 
considered,  the  integral  boundary-layer  method  by  Whitfield  and  the  finite  difference 
boundary-layer  method  of  Herring  and  Mellor  with  algebraic  turbulence  modeling 
provided  the  best  agreement  with  the  nonseparated.  transonic,  adiabatic  wall  data 
considered.  These  methods  were  also  the  most  simple  to  use  and  required  less  computer 
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resources  than  the  other  methods  considered.  This  result,  however,  must  be  tempered  by 
the  fact  that  all  calculation  procedures  did  not  require  and  could  not  be  started  with 
identical  initial  conditions. 
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NOMENCLATURE 

c Airfoil  chord;  also  [3.1  - 1.82  Re^x  x 106]  used  in  Eq.  (25) 

Cf  c Local  skin  friction  coefficient  based  on  boundary-layer  edge  conditions, 

2iWjOeuc2 

cp  Specific  heat 

Cp  Pressure  coefficient  = 2(p  - pJ/p^Uj2 

C£  Critical  pressure  coefficient 

D Dissipation  integral,  Eq.  (5) 

e Turbulent  kinetic  energy  (TKE)  per  unit  mass 

ei  Local  mean  internal  energy 

et  Local  mean  total  energy,  p[ej  +(u2+v2)/2] 

e^.,e^  Unit  vectors  in  the  x and  y directions,  respectively 

h Local  mean  fluid  enthalpy,  e,  + .p/p 

H Local  boundary-layer  shape  factor,  8*/6 

H Column  vector  defined  by  Eq.  (35) 

k Local  mean  thermal  conductivity;  also  exponent  governing  whether  flow  is 

taken  to  be  planar  (k  = 0)  or  axisymmetric  (k  = 1) 

£ Mixing  length,  Eqs.  (42)  and  (43) 

L Length  of  body;  also  dissipation  length  used  in  Eq.  (25) 

M Mach  number 

n Unit  vector  normal  to  a grid  cell  face 

p Local  mean  pressure 

Pr  Prandtl  number 

Prt  Turbulent  Prandtl  number 
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q Local  mean  heat  flux 

r Transverse  wall  radius  for  axisymmetric  flow 

R Local  Reynolds  number  defined  by  Eq.  (17) 

Re  Reynolds  number,  pus  In 

Rs  Longitudinal  radius  of  curvature  of  surface  used  in  Eqs.  (32)  and  (33) 

S Length  variable 

t Time 

T Local  mean  temperature 

u Local  mean  fluid  velocity  component  parallel  to  body  surface 

U Column  vector  defined  in  Eq.  (35) 

Ue  Local  mean  boundary-layer-edge  velocity 

v Local  mean  fluid  velocity  component  normal  to  body  surface 

V Local  mean  velocity  vector,  uex  + vey 

x Axial  or  streamwise  coordinate 

y Coordinate  normal  to  body  surface 

y+  Law-of-the-wall  coordinate, 

a Airfoil  angle  of  attack;  also  variable  used  in  Eq.  (33) 

a*  Variable  used  in  Eq.  (32) 

Oef  Local  mean  effective  thermal  conductivity  defined  by  Eq.  (15) 

0 Angle  between  tangent  to  body  surface  and  free-stream  velocity  vector,  Eq. 

(24);  also  constant  used  in  Eq.  (33) 

0*  Constant  used  in  Eq.  (33) 

S Local  boundary-layer  thickness 
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6*  Local  boundary-layer  displacement  thickness 

5**  Defined  by  Eq.  (3) 

5*  Incompressible  displacement  thickness,  Eq.  (8) 

?T*  Local  boundary-layer  displacement  thickness  which  includes  body  radius 

term,  Eq.  (44) 

e Eddy  viscosity;  also  eddy  diffusivity,  Eqs.  (29  through  33) 

8 Local  boundary-layer  momentum  thickness 

8*  Defined  by  Eq.  (4) 

Local  boundary-layer  momentum  thickness  which  includes  body  radius  term. 
Eq.  (45) 

p Local  mean  molecular  viscosity 

v Local  mean  kinematic  viscosity,  p/p 

I 

yef  Local  mean  effective  kinematic  viscosity  defined  by  Eq,  (14) 

p Local  mean  fluid  density 

a, a*  Constants  used  in  Eq.  (33) 

ot  Schmidt  number  for  turbulent  kinetic  energy 

ax,oy  Defined  by  Eqs.  (38) 

r0  Local  total  shear  stress  (molecular  + turbulent) 

Shear  stress  tensor  defined  by  Eq.  (37) 
rxy,Tyx  Defined  by  Eq.  (40) 

0 Defined  by  Eq.  (18) 

<t>  Defined  by  Eq.  (19) 

X Defined  by  Eq.  (18) 

X Defined  by  Eq.  (16) 
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\}i  Variable  used  in  Eq.  <33) 

qj  Turbulent  dissipation  rate 

SUBSCRIPTS 
aw  Adiabatic  wall 

c Airfoil  chord  length 

e Edge  condition 

local  Local  condition 

m Denotes  middle  region  of  boundary  layer 

o Total  condition;  also  denotes  outer  region  of  boundary  layer 

w Wall  condition 

s Based  on  body  length 

' Undisturbed  free-stream  condition 

6*  , Based  on  boundary-layer  displacement  thickness 

6 Based  on  boundary-layer  momentum  thickness 

SUPERSCRIPTS 

' Denotes  fluctuating  quantity 

— Denotes  nondimensional  quantity,  except  8* 

SPECIAL  NOTATION  AND  CLARIFICATION 

< > Denotes  time  average,  e.g.,  <uv>  = limit  t -*  “ l/2t  jrl(uv)dt'  The  variables 

denoted  in  the  nomenclature  as  mean  quantities  are  understood  to  be  time 
averaged. 
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